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ABSTRACT 

We have obtained broad-band near-infrared photometry for seven Galactic star clusters (M92, M 15, M 13, M5, 
NGC 1851, M 71 and NGC6791) using the WIRCam wide-field imager on the Canada-France-Hawaii Telescope, 
supplemented by images of NGC 1851 taken with HAWK-I on the VLT. In addition, 2MASS observations of the 
[Fe/H] w 0.0 open cluster M67 were added to the cluster database. From the resultant (V -J)-V and (V -Ks)-V 
colour-magnitude diagrams (CMDs), fiducial sequences spanning the range in metallicity, -2.4 < [Fe/H] < +0.3, 
have been defined which extend (for most clusters) from the tip of the red-giant branch (RGB) to '--^ 2.5 magnitudes 
below the main-sequence turnoff. These fiducials provide a valuable set of empirical isochrones for the interpre- 
tation of stellar population data in the 2MASS system. We also compare our newly derived CMDs to Victoria 
isochrones that have been transformed to the observed plane using recent empirical and theoretical colour- Teff 
relations. The models are able to reproduce the entire CMDs of clusters more metal rich than [Fe/H] w -1 .4 quite 
well, on the assumption of the same reddenings and distance moduli that yield good fits of the same isochrones 
to Johnson-Cousins BV{RI)c photometry. However, the predicted giant branches become systematically redder 
than the observed RGBs as the cluster metallicity decreases. Possible explanations for these discrepancies are 
discussed. 



1. INTRODUCTION 

The Two Micron All Sky Survey (2MASS) has uniformly 
scanned the entire sky in three near-infrared bands, JHKs- Cat- 
aloging more than 300 million point sources, 2MASS obser- 
vations are not only advancing our understanding of stars, but 
they also provide the photometric standards used to calibrate all 
other JHKs photometry. In recent years, a new generation of 
JHKs equipped observing facilities have come online, driven 
by the benefits of observing at these longer wavelengths. In 
particular, reduced dust attenuation allows the near-infrared to 
be a better probe of stellar populations in dust-obscured and 
heavily reddened galaxies. Furthermore, the integrated lumi- 
nosities of intermediate-age and old stellar populations, which 
are predominantly due to their asymptotic-giant-branch (AGB) 
and red-giant-branch (RGB) stars, are brightest at these longer 
wavelengths. 

As the infrared (IR) region of the electromagnetic spec- 
trum receives growing attention in modern astrophysics, it be- 
comes desirable to have deep near-IR photometry for Galac- 
tic star clusters. Currently, the available near-IR fiducials of 
these systems remain restricted in both metallicity and lumi- 
nosity, with the observations rarely extend ing to fainter mag- 
nitudes than the base of the RGB (e.g., iFrogel et al.l 119831; 
IValenti. Ferraro. & Ori glia'2007l: lFerraro et al.l200d r^ar clus- 
ters are the ideal stellar populations to observe because, notwith- 
standing a growing number of exceptions (e.g., u Cen), their 
stars are or can be considered homogeneous in both age and 
initial chemical composition to a rather good approximation, at 
least for a number of purposes. Therefore, observations of these 
systems provide us with exceedingly valuable stellar popula- 
tion templates across a broad range of stellar parameter space, 
as well as the data which are needed to test and refine the pre- 
dicted colours from model atmospheres and the temperatures 



given by stellar evolutionary models. 

Stellar population studies rely on the accuracy of both colour- 
Teff relations and the stellar Teft scale in order to derive the 
ages, metallicities, and star formation histories of stellar sys- 
tems based on isochrone fits to observed colour-magnitude di- 
agrams (CMDs). One way to test these relations is to obtain 
photometry in many bandpasses and then to determine the level 
of consistency across all possible CMDs (and colour-colour 
diagrams) that can be generated. While considerable work 
has been carried out to test and improve the colour transfor- 
mations for the BV(RI)c and Stromgren filte r systems (e.g., 
IVandenBerg & Clemj llool IClem et al.1 l2004h . very httle has 
been done to date on the colour- Teff relations f or the near-IR 
(two of the few works to date are Pinsonneault et al. 1 (12004 
and citetGS2003). 

In this paper, new CMDs for open and globular star clus- 
ters, along with published data for field subdwarfs, are com- 
pared with isochrones in order to assess the reliability of re- 
cent (y-y)-reff and {V - Ks)-Tf.ff relations. Additionally, and 
perhaps most importantly, we can for the first time evaluate 
the consistency of the isochrone fits across optical and infrared 
colours. The companion work to this paper, VandenBerg et al. 
(AJ, submitted); hereafter known as VCSIO), focusses on the 
colour-reff relations for the Johnson-Cousins BV(RI)c photo- 
metric system. In both papers, two recently developed sets of 
colour transformations are tested. The first of these has been de- 
rived (see VCSIO) from the latest MARCS model atmospheres 
(Gustafsson et al. 2008), w hile the second is the set of empiri- 
cal relations developed by ICasagrande et al.l (1201 Ol) who used 
the Infrared Flux Method (IRFM) to produce colour-Feff re- 
lations for dwarf and subgiant stars spanning a large range in 
metallicity. 
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Addressing the need for precise fiducial sequences in the 
near-infrared, we have obtained observations of the seven Galac- 
tic stai- clusters, M92, M15, M13, M5, NGC1851, M71 and 
NGC6791. In ^ we describe our observing programs and 
present salient details concerning the reduction of the data — 
including, in particular, the calibration of our photometry to 
the standard 2MASS system. In jTTI the CMDs and clus- 
ter fiducials are presented which provide template stellar pop- 
ulation sequences for the range in metallicity, -2.4 < [Fe/H] 
< +0.3. An examination of how well iso chrones that em- 
ploy th e synthetic MARCS, and the empirical lCasagrande et al.l 
( I2OIOI) . colour- Feff relations are able to reproduce the observed 
VJKs photometry of the local subdwarfs and of our target clus- 
ters is presented in ^ Finally, a short summary of our results, 
as well as a brief discussion of the usefulness of the derived 
fiducial sequences for stellar populations research, is given in 
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Fig. 1 . — Aj-band photometric differences for star.s observed in both 2MASS 
and our VLT fields for NGC 1851. 



2. OBSERVATIONS 

In the 2008B semester, we received WIRCam time on the 
CFHT to observe five Galactic star clusters (NGC 679 1, M 13, 
M 15, M92 and M71) in J and Ks. In the following semester 
(2009 A), we obtained observations of M5, but only in the J 
band due to poor weather conditions. In addition, A'5-band im- 
ages of NGC 1851 that had been previously collected by two 
of us (G.B. and M.D.) using the VLT HAWK-I detector were 
included in the dataset. 

All of our CFHT observations aimed to reach a signal-to- 
noise ratio of 25 at ^ 2.5 magnitudes below the main sequence 
turn-off (MSTO) of each cluster To accomplish this, a series 
of long exposures, together with several short-exposure images 
(which were critical for the calibration of our fields to 2MASS 
photometry), were taken in the J and Ks filters for each cluster 
on our target list. By taking a series of exposures, each star was 
detected multiple times, thereby helping to improve the preci- 
sion of the final photometry. We direct the reader to Table 1 



which lists the image properties of each cluster. Unlike opti- 
cal photometry, where one would separately observe standard 
star fields, the 2MASS All-Sky Point Source Catalog contains 
enough stars in each of our clusters which can be used as stan- 
dai-ds. For example, M 13 has 3221 stars in the 2MASS All-Sky 
Point Source Catalog within 30 arcmin of its center. Of these, 
we selected stars with the lowest claimed photometric errors 
(typically errors < 0.02) for the calibration of our frames. 

In addition to the on-target images, it was necessary to obtain 
a large number of off-target images of the sky because both the 
spatial and temporal variations of atmospheric emission in the 
IR are quite considerable (^ 10% over a timespan of as little as 
10 minutes). Indeed, particular care was taken in determining 
how to construct and subtract a sky image from a target image, 
since the quality of a processed image is dominated by pre- 
cisely how this subtraction is carried out. In order to minimize 
the effects of the variable sky on our photometry, we decided 
to image an equal number of sky frames off-target as science 
frames on-target. We also chose to employ a large dither pat- 
tern so as to remove bad pixels when median combining the 
frames. Once observed, images were pre-processed by CFHT 
staff using the WIRCam pipeline. This included flat fielding, 
bias and dark subtraction, as well as sky subtraction. 

Upon receiving the pre-processed images, instrumental mag- 
nitudes for all stars were obtained using the point spread func- 
tion (PSF) modeling and fitting techniques in the DAOPHOT/ 

ALLS TAR/ ALLFRAME packages ( Stetson 1987, Stetson & Ha^ 
[Ml). In essence, these programs work by detecting stars on 
a specific image, building a model PSF from a few isolated, 
bright stars and then subtracting this PSF from all stars de- 
tected. (For a m ore detailed de scription of how how these pro- 
grams work, see IStetson|[T987i ) The standard star calibration 
directly to 2MASS stars in the same frame, was done by solv- 
ing for the zero point, and colour terms which accounts for the 
difference in the central wavelengths of the CFHT and 2MASS 
7, and K filters. 




Fig. 2. — The derived V-(V -J) fiducial sequence of NGC6791 overlaying 
CFHT (filled circles) and 2MASS (open circles) photometry. 
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Cluster 


Dates of Observations 


Exposure time in J 


Exposure time in K 


Weather Conditions 


M 15 


Aug 10-18, 2008 


AC A r\ C '~l 

45s X 40, 5s X 7 


20s X 97 


photometric 


M 09 
iVl yl. 


Ann in 18 9008 


ZjS X J, ^JS Xi^,jSXZi 


9nc V 1 97 Y 1 Q 
ZWS X iZ/, JS X Vy 


photometric 


M13 


Aug 11-12, 2008 


45s X 25 


20s X 153 


photometric 


M5 


June 11,2009 


20s X 56 




photometric 


NGC 1851 


Jan 30, 2008 




10s X 139 


photometric 


M71 


Aug 10-18,2008 


5s X 22 


5sx45 


photometric 


NGC 6791 


Aug 10-18,2008 


5sx21 




photometric 



Table 1 

Description of observations for each target. 



One can check the overall quality of our transformed magni- 
tudes by comparing stars in common between our fields and the 
2MASS catalog. One example of such a comparison is shown 
in Figure [U for NGC 1851, where the differences between the 
standard 2MASS magnitudes and our final calibrated ones are 
plotted against both magnitude and colour Reassuringly, the 
horizontal lines corresponding to zero difference appear to pass 
through the densest concentration of points in all plots. More- 
over, there seem to be no strong systematic trends as a function 
of colour that would indicate the need for higher order colour 
terms in the photometric solutions. 

To produce the CMDs, our final pho tometry for each clus- 
ter wa s combined with V-band data (see lStets"onll20 00; Stetson 
I2OOI . Our CMD of NGC 1851, which is from VLT observa- 
tions, extends from near the RGB tip to ~ 3 magnitudes below 
the MSTO. For the remaining six clusters, which were observed 
using the CFHT, we were able to obtain observations only for 
magnitudes fainter than the base of the RGB. The upper giant- 
branch stars were saturated in all of our frames, despite observ- 
ing the target clusters with the shortest exposures possible on 
WIRCam, This problem occurred mainly because the seeing 
was 0.3 to 0.5" better than we had requested during the nights 



of observation, which resulted in a more concentrated PSF 

Thus, in order to populate the RGBs of the clusters that 
were observed using the CFHT, we queried the 2MASS cata- 
log for all of the stars within 30' of the center of each cluster. 
With a photometric sensitivity of 10 sigma at 7= 15.8 and Ks= 
14.3 mag, 2MASS observations extend at least to the base of 
the giant branch for each of our target clusters. The 2MASS 
photometry for the se lected giant-branch stars in each system 
(ISkrutskieetal]|2006l) wa s then c ombined with V-ba nd pho- 
tometry for the same stars (IStetsonl l 2000; St etsonll2005[) . 

Although it is ideal to have homogeneous observations for 
the entire range in cluster magnitude, the zero-points of our 
WIRCam photometry were set using 2MASS observations from 
the same catalog that was used to identify the RGB stars. There- 
fore, in principle, one should expect that there are no differ- 
ences between these observations. However, in practice, the 
possibility that uncertainties in the zero-points may be appre- 
ciable should be kept in mind when employing the resultant 
CMDs and fiducial sequences. A list of our observed clusters 
and our adopted reddening and distance moduli are given in 
Table 2. 

2.1. Fiducials 





Fig. 3.— The derived V-(V-Ks) and V-(V-J) fiducial sequences of M71 
overlaying CFHT (filled circles) and 2MASS (open circles) photometry. 



Fig. 4. — The derived V-(V-Ks) fiducial sequence of NGC 1851 overlaying 
HAWK-I Ks photometry. 



4 



Brasseur et al. 



> 




(V-J) 

Fig. 5. — The derived V-(V -J) fiducial sequence of M5 overlaying CFHT 
(filled circles) and 2MASS (open circles) photometry. 




1 2 3 4 1 2 3 



(V-K,) (V-J) 

Fig. 6.— The derived V-(V-Ks) and V-{V-J) fiducial sequences of M13 
overlaying CFHT (filled circles) and 2MASS (open circles) photometry. 



Fiducial sequences are ridge lines of the stellar loci in colour- 
magnitude space. The definition of these fiducials from cluster 
photometry is often based on visual inspection of the CMD, 
since automated scripts typically give poor results in regions of 
low star counts and where the magnitude varies weakly with 
colour, such as the subgiant branch (SGB). Moreover, contam- 
ination from field, AGE, binary, and horizontal-branch (HB) 
stars can significantly skew the the computed locus. For these 
reasons, we have derived all fiducials by dividing the magnitude 
axis into small bins (typically ^ 0.15 mag wide, or smaller in 
regions of nearly constant magnitude) and then calculating the 
median colour of those stars which we judge to belong to the 
cluster 

We present our observed CMDs in Figures|2H8]with 2MASS, 
CFHTAVIRCam, and VLT/HAWK-I photometry plotted as open 
circles, filled circles, and open triangles, respectively. Overlaid 
on the colour-magnitude diagrams are the derived fiducial se- 
quences spanning the MS, SGB, and RGB (see Tables 2-9). In 
the case of M 13, M92, and M 15, there were insufficient points 
to define the fiducial for the subgiant branch; consequently, the 
latest Victoria-Regina Models by VandenBerg et al. (in prepa- 
ration) of the same metallicity of each cluster were used to de- 
fine the transition from the MS to the lower RGB in these cases 
(such points are marked with an asterisk in the tables). 

3. TESTING THE COLOUR- Teff RELATIONS USING OBSERVATIONS 
OF FIELD SUBDWARFS AND STAR CLUSTERS 

Transforming isochrones from the theoretical log Teff-M/jo; 
plane to an observed CMD is accomplished through the use 
of colour- Feff relations to link the fundamental stellar parame- 
ters to photometric indices. In this section, we investigate two 
recently developed colour transformations. The first of these 
is an empirical relation developed by .Casagrande et al., (,2010l) 
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Fig. 7.— The derived V-iV-Kg) and V-{V-J) fiducial sequences of M15 
overlaying CFHT (filled circles) and 2MASS (open circles) photometry. 
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Cluster 


Type 


a 


5 


[Fe/H] 


(m-M)v 


E(B-y) 


M15 


globular 


21:29:58 


+ 12:10:01 


-2.4 


15.29 


0.108 


M92 


globular 


17:17:07 


+43:08:12 


-2.4 


14.62 


0.023 


M13 


globular 


16:41:41 


+36:27:37 


-1.60 


14.40 


0.016 


M5 


globular 


15:18:36 


+02:05:00 


-1.40 


14.45 


0.038 


NGC 1851 


globular 


05:14:06 


- 40:02:50 


-1.40 


15.50 


0.034 


M71 


globular 


19:53:46 


+18:46:42 


-0.80 


13.78 


0.220 


NGC 6791 


open 


19:20:53 


+37:46:30 


+0.3 


13.57 


0.150 



Table 2 

Properties of the Galactic star clusters in our survey. 



(hereafter referred to as CRMBA) who have applied the IRFM 
to a large sample of dwarf and subgiant stars of varying metal- 
licity. They present their results as a set of polynomials which 
relate many photometric indices to Teft and [Fe/H]. This em- 
pirical colour- reff relation has the advantage of being largely 
model independent, but it has the limitation of being applicable 
only to dwarf and SGB stars. 

The second relation was derived by one of us (L.C.) from 
syntheti c spectra base d on the latest MARCS model atmo- 
spheres (iGustafsson et al. 2008) using the procedures described 
in VCSIO, with the 2MASS filter transmission curves, absolute 



calibration and zero-points reported in Casagrande, Portinar i & Flvnn 
( 120061) . The MARCS models consist of plane-parallel (for 
dwarf and SGB stars) and spherical (for giant stars) line-blanketed, 
flux-constant stellar atmospheres for large ranges in effective 
temperature, surface gravity, and chemical composition. By 
convolving the spectra derived from these atmospheres with 
the appropriate filter transmission functions, large tables were 
produced that provide the synthetic magnitudes as functions of 
[Fe/H], Teff, and \ogg. Theoretical relations such as these have 
the advantage over the CRMBA relation that the desired pho- 
tometric indices for any given star or stellar model can be ob- 




FIG. 8.— The derived V-(V-Ks) and V-(V-J) fiducial sequences of M92 
overlaying CFHT (filled circles) and 2MASS (open circles) photometry. 



tained simply by interpolating it in the grid of synthetic colours 
for its [Fe/H], logg, and Teft values. However, their accuracy 
relies heavily on whether or not the synthetic spectra are able 
to reproduce the observed spectra of stars, and on how well 
the filter transmission functions and zero-points are defined. It 
should also be noted that we apply the slight modifications to 
the 2MASS model calibration noted in the recent CRMBA pa- 
per. These very slightly affect the JHKs model colours { J ^ J 
- 0.017; H ^ H + 0.016; Ks -^Ks + 0.003) and indeed these 
improve the fits to the data. 

In the following analysis we transform the latest Victoria 
]sochrones ( VandenBerg et al. in preparation) using the MARCS 
and CRMBA colour- Teff relations and test how well they are 
able to reproduce the observed VJKs colours of field subdwarfs 
and our cluster photometry. Whenever possible, distances de- 
rived from Hipparcos parallaxes, along with the most up-to- 
date estimates of the cluster metallicities, have been assumed. 
However this study is much less concerned with the absolute fit 
of the isochrones to the observed CMDs than with the extent 
to which a consistent interpretation of the same cluster can be 
found on different colour planes. Since any inaccuracies in dis- 
tance moduli or [Fe/H] should be evident in all colour planes, 
they should not effect our overall conclusions (though the pos- 
sibility that chemical abundance anomalies may be present, and 
that they might affect the fluxes in some filter bandpasses more 
than others should be kept in mind). This does, of course, have 
the caveat that the observations must be free from systematic 
errors and that the extinction in all colour bands is accurately 
determined. 

As previously mentioned, this paper provides an extension of 
the BV(7;/)c study by VCSIO to the near-lR. Consequently, we 
have tried to ensure that the analysis of the colour transforma- 
tions for a given cluster is consistent across optical and infrared 
colours by employing identical isochrones (i.e., for the same 
age and metallicity) and by adopting the same values of the 
reddening and distance modulus (in the case of clusters which 
are in common with VCS 1 0: M 67, NGC 679 1 , NGC 1 85 1 , M 5, 
M92). Additionally, the V-band photometry of the clusters pre- 
sented in the following secti ons is from the same catalogues 
( IStetsonll2000l: IStetsonlllOOSh that were used for the compar- 
isons presented by VCSIO. 

The extinction coefficients adopted in this work are E(y - 
y)/E(B-y}=2.25, E{y-K s)IE{B-V)=2n6 and E(V-H)IE{B- 
V)=2.55 dMcCaUi (l2004 ). 

3.1. M67 ([Fe/H] fa 0.0) 

M67 is a particularly well studied, relatively nearby open 
cluster, with a metallicity that is very close to solar, accordi ng to 
the results of high-resolution spectroscopy (.Tautvaisiene et al. I 
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I2OOOI: iRandich et ani2006l) . See VCSIO for a summary of its 
basic properties, together with supporting references. Because 
of its proximity, 2MASS JHKs observations of M67 reach 
well below the MSTO; consequently, the CMDs shown in Fig- 
ure |9]could_be_produ£edsim£ly by combining M67 data from 
2MASS (ISkrutskie et alj|2006h with the latest reduction of the 
V-band photometry discussed bv lStetsonI (l2005h . The isochrone 
in this figure is that for an age of 3.7 Gyr from Micha ud et aTl 
(|2004), and just as VCSIO found from their comparisons of the 
same isochrone with BV(RI)c photometry, both the MARCS 
and CRMBA transformations to the V -J and V -Ks colour 
planes enable the models to match the observed photometry 
exceptionally well over the entire range of luminosity that has 
been plotted. Indeed, the differences between the two colour- 
Teff relations are barely discernible. (Only in the case of M67 
have we compared theory and observations on a [{V-H)q, My]- 
diagram. Although the isochrone fits the data quite well, there 
is a slight, apparently nearly constant offset between the two, 
which may be due to a small error in our assumed value of 
E{V-H) or in the zero-point of the H magnitudes derived from 
the MARCS model atmospheres.) 

3.2. NGC6791 ([Fe/H] ^ +0.3) 

Being the most metal rich cluster in our sample (and only 
open cluster for which we have collected JKs observations), 
NGC6791 provides an important test of the colour- reff re- 
lations applicable to super-metal-rich stars. Here we have 
adopted the recent [Fe/H] estimate from Boesgaard et al. I ( |2009l) . 
This is consistent with iBrogaard et al. I ( 20101) who derived 
[Fe/H] = 0.29 ± 0.03 (rando m) ib 0.07 (sys t emati c) in their 
spectro scopic analysis . Anthon v-Twarog et al.l (l2007l) have noted 
that the lCarraroet al. I (1200 6) determination of [Fe/H] = H-0.38 
for NGC 6791 should be reduced to ^ H-0.28 (hence in very 
good agreement with the Brogaard et al. result) if those authors 



had adopted E(B-V) = 0. 15 in their analysis, instead of 0.09. 

Because of its large distance and low Galactic latitude, red- 
dening estimates of N GC 6791 var y considerably in the litera- 
ture. As discussed bv IChabover. Gr een. &Liebert (1999), the 
derived red denings for NGC 6791 span the range 0.09 < E{B- 
V) < 0.26. iKuchinskietaLl (1 19951) have used subdwarf-B stars 
to provide a tight constraint on the red dening, finding ^(g-y ) = 
0.17 ±0.01. This agrees well with the lSchlegel et all (11998') es- 
timate of E(B-V) = 0.155 mag, which is also favoured by the 

comparison of the NGC 6791 CMD with that for solar neigh- 

borhood stars from Hipparcos data (see'Sand age. Lubin. & VandenBerd 
2003). Additionally, Brogaard et al. (2010) also derived E(B- 
V) = 0.160 ± 0.025, which adds to the support of our adopted 
reddening of 0.15. 

In consistency with VCSIO, we adopt a distance modulus of 
(m-M)v = 13.57 which was found through a fit of the cluster 
CMD to local field dwarfs having metal abundances in the range 
+0.15 < [Fe/H] < +0.45 and criMy) < 0.15 (based on Hippar- 
cos parallaxes). Figure [TOl shows our [(y-7)o,Mv] CMD for 
dwarf and SGB stars in NGC 6791 is well reproduced by ei- 
ther the MARCS- or CRMBA-transformed Victoria isochrone 
on the assumption of the aforementioned reddening and dis- 
tance. (Only along the upper RGB does the isochrone seem to 
deviate to the blue of the observations, which is consistent with 
the possibility that the MARCS atmospheres have insufficient 
blanketing in cool super-metal-rich stars.) Thus, there is excel- 
lent consistency with the isochrone fits reported by VCSIO to 
the [(B-y)o,My]- and [(y-/)o, Mv]-diagrams o f NGC 6791. 
Since E(V-J) = 2.25 E(B-V) (e.g..'M cCall1l2004 . this consis- 
tency provides a further (strong) argument that the foreground 
reddening must, indeed, be quite close to E(B-V) = 0.15. This 
example thus provides an instructive example of how valuable 
it is to have both optical and near-IR photometry to constrain 
the reddening in the case of highly reddened systems. 
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Fig. 9. — Comparison of an 3.7 Gyr isochrone for [Fe/H] = 0.0 with 
2MASS/Stetson V-Ks,V-J, and V-H obs ervations of M67 The so lid and 
dot-dashed curves assume the MARCS and ICasagrande et all 120101) colour 
transformations, respectively. 



Fig. 10. — Comparison of an 8.0 Gyr isochrone for [Fe/H] = +0.3 and 
Y = 0.30 with our V -J observations of NGC6791. The solid and dot-dashed 
curves assume the MARCS and iCasagrande et al., tlOlOi) colour transforma- 
tions, respectively. 
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3.3. M71 ([Fe/H] M-0.75) 

M71 represents the metal -rich end of our sample of globu- 
lar clusters (GCs). We adopt an [Fe/H]=-0.80 as derived by 
iBoe sgaard e t al. I (l2005h . It is quite obvious from the CMD in 
Figure[3]that this cluster has a low Galactic latitude, as contam- 
ination from field stars is still significant even when only those 
stars contained within an annulus of r^l a rcmin around the clus- 
ter center are plotted. Furthermore, the ISchlegel et al.l (Il998l) 
dust maps indicate that the reddening is differential across the 
cluster, which contributes to the spread in colour at any given 
magnitude, thereby hindering the definition of tight photomet- 
ric sequences for this cluster in any photometric band. 

As far as the cluster distance is concerned, we have adopted a 
distance modulus that leads to the most consistent interpretation 
between M71 (in this paper), and that of 47 Tuc (in VCSIO), 
which is known to have very close to the same metallicity. If 
the latter has (m-M)v = 13.40, as assumed by VCSIO, and if 
the diffe rence in magn i tude o f the HBs of 47 Tuc and M7 1 is as 
given bv lHodder eld] ( |1992|) . then a value of {m-M)v = 13.78 
is obtained for M 7 1 . 

Because of its low Galactic latitude (b =-4.6°), M71 poses a 
challenge for reddening determinations: estimates in the liter- 
ature range from E{B-V)= 0.21 to 0.32 (e.g. Kron & Guetter 
1976; Harris 1996). 

Figure [TT] shows the isochrones compared with our obser- 
vations. The observations are matched well by isochrones on 
the [{V-Ks)o, Mv] -plane using either the MARCS or CRMBA 
colour transformations. In the case of the V-J observations, the 
isochrones match the base of the RGB well, but drift to the blue 
side of the lower MS, and to the red of the upper MS. However, 
the discrepancies at the faintest magnitudes may have an ob- 
servational origin given that they arise close to the photometric 
limit of J. It should be noted that the reddening adopted here, 
E{B-V) = 0.20 was chosen to obtain a fit to the B-V photom- 



etry. Thus, the models provide consistent interpretations across 
the B-V,V -J and V -K colour-planes. 

3.4. M5 and NGC 1851 ([Fe/H] « -1.4) 

As discussed quite extensively by VCSIO, when Schlegel et 
al. (1998) reddenings are applied to the observed stars in M5 
and NGC 1851 and then their CMDs are shifted in the verti- 
cal direction so that the red HB populations of the two clusters 
have the same luminosities, one finds that the main-sequence 
fiducials of both clusters superimpose one another nearly per- 
fectly, suggesting a common metallicity. This conflicts with 
the results of most spectroscopic studies, which indicate that 
NGC 1851 i s » 0.2 dex more me tal-rich than M5 (e.g., Carretta 
et al. 2009; iKraft & IvansI (l2003h ). If this is coiTect, the latter 
should have a slightly brighter HB than the former, and there- 
fore a slightly brighter MS at a given color — but this would be 
inconsistent with the assumed difference in [Fe/H] (unless Y or 
the heavy-element mixture also varies). However, a resolution 
of this issue is not important for the present investigation, which 
seeks to determine if consistent interpretations of the photom- 
etry across all colour planes can be found. Following VCSIO, 
we adopt [Fe/H]=-1.4 for both M5 and NGC 1851 which is the 
m etallicity which te nds to be favoured by stellar models such 
as lVandenBergl (2000) (see VCSIO for a detailed description on 
this choice of metallicity). 

In Figure [12] we overlay the same 1 1 Gyr isochrone used 
by VCSIO in their study of BV(RI)c photometry, onto the 
[(y-/rs)o, My] -diagram of NGC 1851. For consistency with 
VCSIO, we have adopted the same re ddening and distance 
modulus; specifically, E(B-V) = 0.034 (ISchlegel et al.lll998l) 
and (m-M )v = 15.50, based on obs erved and predicted HB lu- 
minosities (VandenBerg et al.ll2000l) as well as MS-fits of clus- 
ter fiducials to the CRMBA-transformed isochrones (and hence 
to the subdwarf standards for which these relations are defined; 



M71 

(m-M)y= 13.78 
E(B-V) = 0.20 
[Fe/H] = -0.80 





NGC1851 
(m-M)^,= 15.50 
E(B-V) = 0.034 
[Fe/H] = -1.4 




Fig. 1 1 . — Comparison of an 11 Gyr isochrone for [Fe/H] = -0.75, with our 
obs ervations of M 7 1 . The solid and dot-dashed curves assume the MARCS 
and lCasagrande et aljj201(i) colour transformations, respectively. 



Fig. 12. — Comparison of an 11 Gyr isochrone for [Fe/H] = —1.4, [a/Fe] 
» 0.4, and y = 0.25 with our V-Ks observat ions of NGC 1851 The so hd and 
dot-dashed curves assume the MARCS and ICasagrande et alj 120101) colour 
transformations, respectively. 
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VCSIO). We find that isochrones using the MARCS colour-7;ff 
relations, or the CRMBA transformations for dwarf/SGB stars, 
are able to reproduce our V -Ks,V photometry of NGC 1851 
rather well, except for the upper RGB (where there are few stars 
and where the stellar images are approaching the saturation 
limit). These results are fully consistent with those found by 
VCSIO when they analyzed V-I,V data. However, as they re- 
ported, the giant branch of the isochrone on the [{B-V){), My]- 
diagram is significantly bluer than the observed RGB. Thus, 
only the MS stars can be consistently fitted on all of the colour 
planes that have been considered. 

In the case of M5, we have assumed that E(B-V) = 0.038 
and (m-M)v = 14.45 (see VCSIO), which results in the com- 
parison of the MARCS- and CRMBA-transformed isochrones 
with the V -J observations shown in Figure [13] (Note that the 
isochrones are identical to those fitted to NGC1851 photome- 
try.) Given the differences in the relative positions of the RGBs 
of M5 and NGC 1851 noted above, it is not surprising to find 
that, since the MARCS models are able to reproduce the giant 
branch of NGC 1851 in V-Ks and V-I, they lie to the red of the 
observed RGB of M 5 in V -7 and V- / (see VCSIO, who also 
show that the predicted giant branch matches the B-V colours 
of M5 giants quite well). As far as the main sequence of of M5 
on the [(y-y)o, Mv] -diagram is concerned, the isochrones lie 
slightly blueward of the lower MS, and to the red of the upper 
MS, as seen in M71. As discussed below, this is also seen in 
the case of M 15 and M92, and perhaps indicates that there is a 
small problem with the transformations toV -J. 

3.5. M13 ([Fe/H]K. -1.6) 

For our comparisons with M 13, we ha ve selected isochrone s 
with [Fe/H] = -1 .6 based on the work of iKraft & ly aiis (2003^ 
This is within 0.05 dex of the long favoured |Zinn & Wes3 



1£84) value of-1 .65 and within 0.02 dex of the recent lCarretta et 
(I2009i) value of -1.5 8. We have also adopt ed (m-M)v = 14.40 



which is very close to the value obtained bv lGrundahl Stetson. & Andersen 



(119981) from a fit of Stromgren photometry to metal-poor subd- 
warfs with Hipparc os parallaxes, who also assumed [Fe/H] = 
-1.6. Additionally, iPaltrinieri et al. I (Il998l) derived (5(V) be- 
tween M3 and M13 of 0.64 mag - so if M3 has 15.00 (as 
adopted in VCSIO), then 14.40 is impHed for M13 from the 
Paltrinieri et al result (if the 0.01 mag difference in redden- 
ing is taken into ac count). We adopt E{B-V) = 0.016 from 
ISchlegeleTaL (1998). 

In figure[T4l we compare a 12 Gyr isochrone with our VJ and 
VKs photometry, finding that (for both colours) the isochrones 
provide a good match to the observed MS, but not to the ob- 
served RGB, which is significantly bluer at a fixed magnitude 
than the predicted colours. This is also found when fitting 
VI data (not shown here); i.e., when the same isochrones are 
overlaid onto V -I,V photometry, the main sequence is well 
matched but the predicted giant branch is offset to the red of 
the observed RGB. Interestingly, no such problem is found fit- 
ting the VI (or BV) photometry of M3 (see VCSIO), which is 
thought to have close to the same metal abundance as M 13. 

3.6. M15 and M92 ([Fe/H] K, -2.4) 

M 15 and M92 constitute the most metal poor clusters in our 
sample. For both clusters, a metallicity near [Fe/H] = -2.4 
is consiste nt with the latest results from high-reso lution spec- 
troscopy (.Kraft & Ivansll2003l:ICarretta et al.ll2009h . As regards 
thei r distances, we have ado pted (m-M)v = 14.62 for M92 
(see|\MdenBerget al.' ('2002')), and {m-M)v = 15.29 for M 15 
given that Durrell & Harris (1993) have found that the main 
sequences of both clusters will superimpose one another if a 
vertical shift of 0.67 mag is applied to the M 92 CMP once 
their t urnoff colours are matched. According to lSchlegel et al.l 
(119981) . the reddenings of M92 and M 15 are E{B-V) = 0.023 
and 0.108 mag, respectively. 
ID As shown by VCSIO, a MARCS-transformed isochrone for 
13.5 G yr and [Fe/H] = -2.4 provides a reasonably good fit to 
D 



(m-M)y= 14.45 
E(B-V) = 0,038 
[Fe/H] = -°1.40 




4 - 



(V-J)o 

Fig. 13. — As in the previous figure, except that the (same) isochrone is 
compared with our V -J observations of M5. 



M13 

(m-M)^,= 14.40 
E(B-V) = 0.016 
[Fe/H] = -1.6 



(V-K,)„ 




(V-J)o 



Fig. 14. — Comparison of a 12 Gyr isochrone for [Fe/H] = -1.60 with our 
observations of M 13. The solid and dot-dashed curves assume the MARCS 
and Casagrande et all 1201 J) colour transformations, respectively. 
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both the BV and VI photometry of M92 — aside from the fact 
that a small colour shift must be applied to the V -I colours 
in order to obtain consistent fits to the various CMDs that can 
be generated from BV(RI)c photometry. Their additional ob- 
servation that the isochrones tend to deviate to the blue side 
of the observed lower MS is apparently seen in especially the 
[(y -7)o, Mi/]-diagrams of M15 and M92 as well (see Fig- 
ures [15] and [161). We have no explanation for such deviations, 
which appear to be most pronounced in the VJ photometry for 
M15. 

Perhaps the most serious concern is that, in both the V -J 
and V -Ks colour planes, the isochrones lie consistently to the 
red of the observed RGBs of M 15 and M92. Near the base of 
the giant branch, the offsets amount to ~ 0.12 mag in V -Ks 
and ^ 0.09 mag inV-J and the typical error in the photometry 
0.02 mag) cannot explain these colour offsets. Surprisingly, 
VCSIO find that the same isochrones provide a good fit to op- 
tical CMDs, apart from the small, constant offset between the 
predicted and observed RGB on the [(V -/)o, My] -plane pre- 
viously mentioned. Since we are using the same V-band pho- 
tometry as VCSIO, one might question the accuracy of the J 
and Ks photometry for M15 and M92. However, this is not 
a viable explanation because the near-IR photometry for the 
cluster giants, as with all our CFHT clusters, comes directly 
from 2MASS, and such a large zero-point offset is highly un- 
likely. Indeed, because the isochrones reproduce BV(RI)c ob- 
servations quite well, one wonders whether there is a problem 
with the synthetic J and Ks magnitudes that are derived from 
MARCS model atmospheres, though we have no other reason 
to question to these results. Further work is clearly needed to 
understand this puzzle. 

3.7. Field Subdwarfs (-2.2 < [Fe/H] < - 0.5) 

Subdwarfs are metal-poor stars with halo kinematics whose 
orbits have brought them into the solar neighbourhood. With 



many of these stars having accurate distances determinations 
from Hipparcos and recent estimates of log g, and [Fe/H] based 
on high-resolution spectroscopy, they provide strong constraints 
on the Teff scale of isochrones and colour- Teff relations at low 
metallicities. In the following analysis, we use a sample of 
^100 subdwarfs compiled by VCSIO. The photometric and 
fundamental stellar properties (e.g., Teff, logg and [Fe/H]) for 
the majority of these stars were taken from Table 8 of CRMBA 
with the addition of 4 stars from .Clementini et al. (1999). For a 
further description on the properties and choice of this sample, 
see VCSIO. 

In figure [17] we show the level of agreement between the ob- 
served and predicted 2MASS colours, (V -Ks), {V -J), and 
(J-Ks), of the subdwarfs in our sample. The predicted MARCS 
colours for each star are found by interpolating in tables of syn- 
thetic colours assuming the T^ff, log g, and [Fe/H] values given 
by CRMBA. For these colour-colour comparisons we use only 
MARCS transformations since the CRMBA relations are them- 
selves based on most of the same stars used here, so they will 
necessarily yield colours that agree well, in the mean, with the 
observed photometry. As noted in the three panels of figure [TT] 
offsets of ^0.01 mag to the blue for the predicted (V-Ks) and 
(J-Ks) are needed to achieve consistency with observations. 
The predicted (V - J) colours show a very small mean offset of 
0.001 mag from the observed colours, though the largest devia- 
tions occur for the reddest stars, which suggests (perhaps) that 
the transformations to V -J for cooler dwarf stars may need 
some adjustment. (Recall that, although this evidence is based 
on very few cool subdwarfs with measured V -J colours, simi- 
lar problems were seen in our comparisons with the V -J pho- 
tometry for MS stars in the low metallicity clusters M92, M 15, 
M 13 and M5; see the previous subsections.) 

While the comparisons described above are completely inde- 
pendent of stellar evolutionary models, instructive comparisons 
of the theoretical isochrones with the Hipparcos subdwarfs can 
also be carried out. By superimposing subdwarfs onto a grid of 




Fig. 15. — Compai'ison of a 13.5 Gyr isochrone for [Fe/H] = —2.40 with our 
observations of Ml 5 The solid and dot-dashed curves assume the MARCS FiG- 16.— As in the previous figure, except that the (same) isochrone is 

and .Casagrande'etai] jMcj) colour transformations, respectively. compared with our observations of M92. 
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isochrones for a wide range in [Fe/H] (and a fixed age — though 
isochrones exhibit essentially no age dependence at faint lumi- 
nosities), it is possible to examine the consistency of the ob- 
served [Fe/H] values and colours with those inferred from the 
isochrones, given the observed temperatures. If a star lies on the 
same isochrone (say, one with [Fe/H]= -1.5) in all of the colour 
planes, one can conclude that the colour- Teff transformations are 
consistent for those colours. If any obvious discrepancies exist 
for different filter combinations, then this is an indication of a 
problem with the colour-temperature relations for one or more 
colour indices (i.e., they do not adequately represent the spec- 
trum of a star in the wavelength ranges spanned by the relevant 
filters). 

In the lower panel of fi gure [TSl we superimpose Hippar- 
cos dvan Leeuwen I ( l2007h ) subdwarfs onto 12 Gyr Victoria 
isochrones for [Fe/H] values from -2.4 to -0.6 (in the di- 
rection from left to right) transformed into the [{V-J)Mv]- 
, [{V-Ks,Mv]- and [(7 -/Ts), My] -planes using the MARCS 
transformations. For each subdwarf, the CRMBA estimate of 
its [Fe/H] value minus the [Fe/H] value of the isochrone on 
which the subdwarf is located is plotted in the middle panel, 
while the upper panel plots the difference between the observed 
subdwarf colour and that of the isochrone which has the same 
metallicity as the star (at the same My). These panels provide 
a measure of the quality of the fit to the observations. For 
example, if the isochrones were too red, most of the points 
would have positive (5[Fe/H] values and negative (^(colour) val- 
ues. Figures similar to [18] and [191 but for the (B-V), (V-R), 
and (V-I) colour planes, are given in the study by VCSIO, 
where the individual subdwarfs are identified. 

The subdwarfs in the V -Ks and V -J colour planes show 
significant relative shifts with respect to the isochrones (i.e.. 



the subdwarfs do not lie on the same position relative to the 
same isochrone in both colours). However, the scatter on the 
(5[Fe/H] and (5(colour) planes is quite symmetric overall in both 
V -Ks and V -J, with most of the redder stars having nega- 
tive (5[Fe/H] and positive (^(colour) values, whereas most of the 
bluer stars have positive (5[Fe/H] and negative (5(colour) values. 
Thus, the isochrones are consistent with the observations, in the 
mean, with only a slight trend in the sense that the isochrones 
are somewhat too blue for the redder stars and too red for the 
bluer stars. 

A different conclusion is drawn regarding the J -Ks trans- 
formations shown in the right panel of Figure [TU the MARCS - 
transformed isochrones imply implausibly low values of [Fe/H], 
which are contrary to the implications from the other colour 
planes. When looking at the (5[Fe/H] and (5(colour) plots in 
the upper panels, one finds that most of the points have posi- 
tive (5[Fe/H] and negative (^(colour) values, indicating that the 
isochrones are too red, in the mean, relative to the observations. 
The J -Ks panels also indicate that there is no obvious trend 
with colour or with metallicity. 

Figure [19] plots the same isochrones as in Figure [18] but 
transformed to the various colour planes using the CRMBA 
colour-reff relations. When comparing the subdwarf colours 
and metallicities with those of the isochrones, we find there is 
no obvious trend with either colour or metallicity, which indi- 
cates that the CRMBA-transformed isochrones agree well, in 
the mean, with the observations, while the absence of trends 
provides a verification that the functional form used in CRMBA 
transformations well represents the (bulk of) stars upon which 
it is built. 

4. SUMMARY 
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.8 1.2 1.6 2 1.2 1.6 2 2.4 2.8 .4 .6 .8 

(V-J)o (V-Kjo (J-K3)o 

Fig. 18. — Lower Panels: Sample of Hipparcos subdwarfs overlying 12 Gyr Victoria isochrones in the [(y-7)o, [(V-Ks)i), My]- and [(J-Ks)o, Mv]-planes. 
Tlie isoclirones plotted here have [Fe/H] values of -2.4 to -0.6 in increments of 0.2 dex, from left to right, and MARCS transformations have been used to transpose 
the models to the various CMDs. To demonstrate the age insensitivity for magnitudes My > 4.5, dashed lines show 10 Gyr isochrones for [Fe/H] = -0.6, -1.4, 
and -2.4. The photometry and metaUicity of each subdwarf is taken from CRMBA, and the error bars represent la uncertainties in the My values as derived from 
Hipparcos parallaxes. Subdwarfs having [Fe/H]< -1.2, as given by CRMBA, are plotted as filled circles, and those with higher metaUicities are plotted as open 
circles. Middle Panels: Plotted as a function of colour, the difference between the CRMBA estimate of [Fe/H] for each star and that inferred from the interpolated 
(or extrapolated when necessary) isochrone that matches its location on the colour-My diagram in the lower panel. The arrows and the niunbers, together with 
the standard deviation (in parentheses), indicate the mean 5[Fe/H] for all of the subdwarfs that have been considered.. Upper Panels: The difference in colour that 
would need to be applied to each subdwarf in order to achieve perfect consistency of its position relative to the isochrones in the lower panel, n this case, the arrows 
and numbers indicate the mean values of <5(colour) for the subdwarf sample. 




.8 1.2 1.6 2 1.2 1.6 2 2.4 2.8 .4 .6 .8 

(V-J)o (V-Kjo (J-K,)o 

Fig. 19. — Similar to previous figure, except that the isochrones are transformed using the colour-Teff relations developed by CRMBA. 
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We have derived fiducial sequences for the Galactic star clus- 
ters M15, M92, M13, M5, NGC1851, M71, and NGC 6791 
(encompassing -2.4 < [Fe/H] < +0.3) using near-lR observa- 
tions obtained with the WlRCam imager on the Canada-France- 
Hawaii Telescope and the HAWK-1 detector on the VLT. These 
fiducial sequences, which are presented in Figure |20] on the 
[(V -J)o,Mv]- and [(y-/rs)o, Mi/]-planes, illusti-ate how the 
photometric properties of these old stellar populations vary as 
a function of [Fe/H]. With spectroscopic metallicity determi- 
nations and relative age estimates that are accurate to within 
±0.25 dex and ±1.5-2 Gyr, respectively, these fiducials can, 
in principle, be used to photometrically determine the age and 
metallicity of resolved stellar systems. Unlike isochrone analy- 
ses, metallicity determinations made through comparisons with 
fiducials are independent of stellar evolutionary models. Hence, 
these fiducials provide a set of empirical isochrones that can 
serve as valuable tools for future stellar population investiga- 
tions involving the 2MASS filters. 

Based on these data, a summary of our results is as follows: 

(1) The predicted V-J, V-Ks, and J-Ks colours given 
by the MARCS transformations for local subdwarfs agree well 
with those observed (i.e., to within ^ 0.01 mag) when the tem- 
peratures and metallicities of the latter are as given by CRMBA. 

(2) While isochrones appear to be too blue when compared 
with lower main sequence stars in metal-poor GCs on the 
[(V-J)o, Ml/] -plane, no significant problems are found match- 
ing their V -Ks colours. 

(3) Our analysis of cluster RGB stars indicates that the isochrones 
and the colour transformations faithfully reproduce the proper- 
ties of metal-rich giants, but not those of lower metallicity. In 
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Fig. 20. — Our cluster fiducials mapped to the [(V - Ks)o, My]- and 
[(V— 7)o, Mv]-planes using tlie reddenings and distance moduli given in Table 
[2] From left to right, corresponding to increasing metalhcity, are M 15 (solid 
black), M92 (dotted black), M 13 (grey short-dashed), NGC 1851 (black long- 
dashed), M5 (black solid line), M71 (black dot-dashed) and NGC 6791 (grey 
long-short dashed). 



fact, the RGB segments of the isochrones become systemati- 
cally redder than the observed RGBs with decreasing [Fe/H]. It 
seems unlikely that errors in the model Teff scale are sufficient 
to cause this problem, because the discrepancies are primarily 
on the V-J and V -Ks colour planes. The same MARCS- 
transformed isochrones are able to reproduce the CMDs of 
M92 derived from BV{RI)c photometry quite well. This seems 
to suggest that there are problems with the predicted / and Ks 
magnitudes at low metallicities, but we are not able to provide 
an explanation for the origin of such difficulties. 
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V 


V-J 


V 


V-J 


19.817 


2.062 


17.411 


1.762 


19.659 


2.001 


17.437 


1.831 


19.515 


1.944 


17.476 


1.897 


19.363 


1.889 


17.485 


1.953 


19.215 


1.840 


17.459 


2.011 


19.067 


1.787 


17.376 


2.043 


18.924 


1.738 


17.263 


2.058 


18.771 


1.692 


17.109 


2.071 


18.633 


1.656 


16.959 


2.091 


18.486 


1.625 


16.809 


2.114 


18.331 


1.595 


16.659 


2.130 


18.172 


1.567 


16.509 


2.152 


18.011 


1.548 


16.359 


2.164 


17.845 


1.529 


16.209 


2.186 


17.676 


1.528 


16.059 


2.211 


17.516 


1.548 


15.909 


2.236 


17.388 


1.588 


15.759 


2.264 


17.357 


1.639 


15.609 


2.299 


17.367 


1.699 







Table 3 

Fiducial sequence for the open cluster NGC 679 1 as shown in Figure[2] 



V 


V-J 


V-Ks 


V 


V-J 


V-Ks 


19.505 


1.805 


2.408 


15.854 


2.080 


2.734 


19.338 


1.757 


2.338 


15.627 


2.120 


2.766 


19.168 


1.712 


2.273 


15.426 


2.152 


2.798 


18.991 


1.670 


2.214 


15.217 


2.177 


2.835 


18.807 


1.635 


2.163 


14.986 


2.209 


2.880 


18.622 


1.604 


2.138 


14.736 


2.247 


2.935 


18.429 


1.576 


2.099 


14.538 


2.274 


2.974 


18.232 


1.555 


2.068 


14.356 


2.306 


3.020 


18.029 


1.542 


2.051 


14.177 


2.342 


3.071 


17.824 


1.542 


2.040 


13.993 


2.381 


3.129 


17.634 


1.563 


2.110 


13.818 


2.422 


3.187 


17.500 


1.611 


2.199 


13.646 


2.465 


3.250 


17.445 


1.673 


2.278 


13.459 


2.516 


3.324 


17.439 


1.745 


2.353 


13.305 


2.562 


3.390 


17.441 


1.814 


2.424 


13.139 


2.648 


3.481 


17.396 


1.882 


2.492 


12.906 


2.733* 


3.619* 


17.267 


1.931 


2.564 


12.693 


2.833* 


3.800* 


17.084 


1.959 


2.606 


12.523 


2.962* 


3.970* 


16.885 


1.977 


2.631 


12.347 


3.137* 


4.164* 


16.688 


1.989 


2.649 


12.212 


3.300* 


4.350* 


16.485 


2.021 


2.667 


12.116 


3.469* 


4.552* 


16.268 


2.035 


2.687 


12.035 


3.717* 


4.883* 


16.060 


2.070 


2.709 


11.997 


4.109* 


5.401* 



Table 4 

Fiducial sequences for the globular cluster M71 as shown in Figure[3J Due to the sparseness of the RGB, those 
points with V <13 are marked by asterisks to indicate they have large uncertainties associated with them. 
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y 




y 




21.95 


2.01 


17.74 


1.42 


21.84 


1.96 


17.54 


1.46 


21.60 


1.84 


17.32 


1.47 


21.37 


1.74 


17.08 


1.48 


21.16 


1.65 


16.83 


1.50 


20.95 


1.56 


16.57 


1.51 


20.75 


1.50 


16.33 


1.54 


20.55 


1.43 


16.05 


1.57 


20.36 


1.36 


15.81 


1.60 


20.16 


1.30 


15.55 


1.63 


19.96 


1.26 


15.27 


1.67 


19.75 


1.21 


15.04 


1.70 


19.53 


1.17 


14.80 


1.74 


19.31 


1.12 


14.57 


1.79 


19.07 


1.08 


14.35 


1.84 


18.83 


1.05 


14.14 


1.89 


18.59 


1.03 


13.91 


1.95 


18.34 


1.05 


13.71 


2.00 


18.14 


1.09 


13.49 


2.06 


18.00 


1.15 


13.29 


2.12 


17.92 


1.25 


13.09 


2.18 


17.89 


1.28 


12.86 


2.27 


17.87 


1.33 


12.68 


2.36 


17.84 


1.37 


12.50 


2.47 


17.80 


1.39 


12.31 


2.59 



Table 5 

Fiducial sequence for the globular cluster M 5 as shown in Figure[5J 



y 


V-Ks 


y 


V-Ks 


21.993 


2.087 


18.806 


1.914 


21.856 


2.020 


18.695 


1.963 


21.729 


1.959 


18.557 


1.998 


21.609 


1.904 


18.401 


2.024 


21.473 


1.844 


18.245 


2.064 


21.344 


1.789 


18.083 


2.081 


21.212 


1.741 


17.917 


2.099 


21.077 


1.696 


17.740 


2.117 


20.938 


1.652 


17.570 


2.136 


20.797 


1.611 


17.410 


2.155 


20.655 


1.573 


17.237 


2.177 


20.510 


1.536 


17.050 


2.202 


20.360 


1.502 


16.885 


2.237 


20.212 


1.472 


16.720 


2.262 


20.052 


1.444 


16.540 


2.292 


19.893 


1.419 


16.342 


2.348 


19.731 


1.402 


16.180 


2.389 


19.564 


1.393 


16.027 


2.420 


19.399 


1.400 


15.857 


2.475 


19.159 


1.457 


15.702 


2.519 


19.063 


1.511 


15.558 


2.568 


18.996 


1.575 


15.403 


2.617 


18.952 


1.644 


15.262 


2.665 


18.923 


1.725 


15.113 


2.708 


18.906 


1.797 


14.964 


2.768 


18.876 


1.868 


14.815 


2.818 



Table 6 

Fiducial sequence for the globular cluster NGC 1851 as shown in Figure|4] 
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y 


V-J 


y-Ks 


y 


V-J 


y-Ks 


22.039 


2.011 


2.606 


17.548 


1.331* 


1.782* 


21.884 


1.939 


2.513 


17.420 


1.360* 


1.814* 


21.740 


1.883 


2.416 


17.271 


1.381* 


1.844* 


21.606 


1.832 


2.344 


17.112 


1.397 


1.868 


21.460 


1.766 


2.267 


16.952 


1.411 


1.888 


21.306 


1.709 


2.176 


16.786 


1.424 


1.907 


21.173 


1.653 


2.120 


16.433 


1.451 


1.945 


21.030 


1.596 


2.042 


16.104 


1.478 


1.974 


20.891 


1.541 


1.988 


15.747 


1.511 


2.031 


20.761 


1.501 


1.931 


15.585 


1.527 


2.050 


20.640 


1.455 


1.868 


15.408 


1.546 


2.072 


20.502 


1.414 


1.801 


15.214 


1.569 


2.104 


20.372 


1.368 


1.760 


15.036 


1.601 


2.145 


20.249 


1.334 


1.712 


14.880 


1.620 


2.173 


20.112 


1.299 


1.662 


14.715 


1.643 


2.215 


19.979 


1.260 


1.620 


14.546 


1.667 


2.249 


19.839 


1.220 


1.578 


14.393 


1.697 


2.298 


19.700 


1.182 


1.518 


14.238 


1.722 


2.342 


19.560 


1.146 


1.491 


14.088 


1.748 


2.379 


19.411 


1.119 


1.454 


13.945 


1.774 


2.415 


19.263 


1.095 


1.419 


13.800 


1.802 


2.475 


19.113 


1.063 


1.397 


13.646 


1.843 


2.529 


18.962 


1.042 


1.369 


13.499 


1.875 


2.573 


18.800 


1.024 


1.342 


13.360 


1.907 


2.628 


18.637 


1.008 


1.320 


13.211 


1.944 


2.689 


18.472 


0.999 


1.327 


13.051 


1.999 


2.768 


18.306 


0.999 


1.338 


12.923 


2.053 


2.818 


18.152 


1.014 


1.368 


12.787 


2.104 


2.880 


18.022 


1.043 


1.430 


12.642 


2.182 


2.952 


17.923 


1.083* 


1.506* 


12.487 


2.258 


3.058 


17.842 


1.139* 


1.556* 


12.289 


2.373 


3.223 


17.751 


1.192* 


1.632* 


12.152 


2.474 


3.352 


17.692 


1.242* 


1.684* 


12.005 


2.612 


3.547 


17.642 


1.290* 


1.734* 


11.895 


2.729 


3.741 



Table 7 

Fiducial sequences for the globular cluster M 13 as shown in Figure^ Due to the lack of sufficient data at some 
magnitudes, isochrones were used to determine the fiducial points indicated by asterisks. therefore these points 

have large uncertainties associated with them. 
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y 


V-J 


y-Ks 


y 


V-J 


y-Ks 


21.994 


1.833 


2.309 


17.541 


1.326* 


1.758* 


21.836 


1.752 


2.241 


17.418 


1.346* 


1.800* 


21.688 


1.684 


2.177 


17.279 


1.368* 


1.832* 


21.550 


1.621 


2.119 


17.130 


1.376* 


1.857* 


21.420 


1.570 


2.064 


16.970 


1.391* 


1.880* 


21.187 


1.474 


1.969 


16.801 


1.406* 


1.901* 


21.064 


1.433 


1.920 


16.623 


1.421* 


1.922* 


20.924 


1.384 


1.867 


16.453 


1.435* 


1.942* 


20.792 


1.347 


1.818 


16.292 


1.448 


1.961 


20.667 


1.314 


1.772 


16.118 


1.463 


1.982 


20.537 


1.280 


1.725 


15.929 


1.479 


2.005 


20.403 


1.247 


1.678 


15.637 


1.507 


2.054 


20.273 


1.216 


1.635 


15.456 


1.525 


2.080 


20.142 


1.186 


1.593 


15.256 


1.547 


2.109 


20.000 


1.155 


1.549 


15.089 


1.565 


2.136 


19.859 


1.127 


1.509 


14.925 


1.595 


2.162 


19.716 


1.099 


1.470 


14.744 


1.627 


2.201 


19.572 


1.073 


1.432 


14.541 


1.653 


2.247 


19.422 


1.047 


1.396 


14.321 


1.692 


2.308 


19.267 


1.023 


1.361 


14.140 


1.727 


2.357 


19.105 


1.000 


1.328 


13.984 


1.761 


2.399 


18.938 


0.980 


1.300 


13.827 


1.783 


2.433 


18.780 


0.965 


1.278 


13.682 


1.818 


2.467 


18.617 


0.953 


1.261 


13.528 


1.857 


2.505 


18.451 


0.951 


1.259 


13.379 


1.889 


2.544 


18.294 


0.963 


1.276 


13.236 


1.920 


2.583 


18.164 


0.990* 


1.314* 


13.085 


1.964 


2.626 


18.062 


1.026* 


1.365* 


12.933 


1.999 


2.675 


17.981 


1.068* 


1.415* 


12.798 


2.046 


2.730 


17.912 


1.113* 


1.450* 


12.655 


2.086 


2.787 


17.852 


1.163* 


1.521* 


12.503 


2.147 


2.854 


17.803 


1.210* 


1.589* 


12.362 


2.208 


2.922 


17.751 


1.261* 


1.662* 


12.235 


2.262 


2.998 


17.652 


1.300* 


1.723* 


12.103 


2.332 


3.108 



Table 8 

Fiducial sequences for the globular cluster M 92 as shown in Figure[8] Due to the lack of sufficient data at some 
magnitudes, isochrones were used to determine the fiducial points indicated by asterisks. therefore these points 

have large uncertainties associated with them. 
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V 


V-J 


V-Ks 


V 


V-J 


V-Ks 


22.010 


1.764 


2.140 


17.979 


1.468* 


1.932* 


21.887 


1.694 


2.079 


17.830 


1.476* 


1.957* 


21.764 


1.653 


2.030 


17.670 


1.491* 


1.980* 


21.624 


1.604 


1.977 


17.501 


1.506* 


2.001* 


21.492 


1.547 


1.928 


17.323 


1.521 


2.022 


21.367 


1.504 


1.882 


17.153 


1.535 


2.072 


21.237 


1.450 


1.835 


16.992 


1.548 


2.081 


21.103 


1.417 


1.788 


16.818 


1.563 


2.122 


20.973 


1.366 


1.745 


16.629 


1.579 


2.145 


20.842 


1.326 


1.703 


16.337 


1.607 


2.204 


20.700 


1.277 


1.659 


16.156 


1.625 


2.230 


20.559 


1.247 


1.619 


15.956 


1.647 


2.259 


20.416 


1.199 


1.580 


15.789 


1.665 


2.296 


20.272 


1.173 


1.542 


15.625 


1.695 


2.302 


20.122 


1.147 


1.506 


15.444 


1.727 


2.331 


19.967 


1.123 


1.471 


15.241 


1.753 


2.367 


19.805 


1.100 


1.438 


15.021 


1.792 


2.418 


19.638 


1.080 


1.410 


14.840 


1.827 


2.457 


19.480 


1.065 


1.388 


14.684 


1.861 


2.499 


19.317 


1.053 


1.371 


14.527 


1.883 


2.523 


19.151 


1.051 


1.369 


14.382 


1.918 


2.557 


18.994 


1.063 


1.386 


14.228 


1.957 


2.595 


18.864 


1.090 


1.424 


14.079 


1.989 


2.644 


18.762 


1.126 


1.475 


13.936 


2.020 


2.683 


18.681 


1.168* 


1.535 


13.785 


2.064 


2.726 


18.612 


1.213* 


1.600 


13.633 


2.099 


2.785 


18.552 


1.263* 


1.671 


13.498 


2.146 


2.850 


18.503 


1.310* 


1.739 


13.355 


2.186 


2.907 


18.451 


1.361* 


1.772 


13.203 


2.227 


2.976 


18.352 


1.400* 


1.833 


13.062 


2.278 


3.062 


18.241 


1.426* 


1.858* 


12.935 


2.322 


3.128 


18.118 


1.446* 


1.900* 


12.803 


2.382 


3.208 



Table 9 

Fiducial sequences for the globular cluster Ml 5 as shown in Figure[7] Due to the lack of sufficient data at some 

MAGNITUDES, ISOCHRONES WERE USED TO DETERMINE THE FIDUCIAL POINTS INDICATED BY ASTERISKS. THEREFORE THESE POINTS 

HAVE LARGE UNCERTAINTIES ASSOCIATED WITH THEM. 



